11 Search Results

This paper describes new modeling capabilities for predicting H-mode pedestal profiles in spherical tokamaks. Temperature profiles for NSTX discharges 132543 and 132588 are modeled by coupling the \textsc{astra} transport solver with neoclassical transport and gyrokinetic-based reduced models for electron temperature gradient (ETG) and kinetic ballooning mode (KBM) instabilities. A quasi-linear surrogate model for ion-scale transport is developed using linear \textsc{gene} simulations, requiring only a single free parameter calibrated to one discharge. Time-evolving the temperatures with fixed density yields good agreement with experiments for both discharges. Systematic analysis of the transport mechanisms reveals that neoclassical transport is huge across the entiremore » pedestal region for the ion channel. ETG turbulence is large in the plasma edge and low density gradient region, contributing substantially to the electron channel. However, KBM/MHD-like modes also drive significant transport in both the ion and electron thermal channels, making them essential for accurate pedestal modeling. Further refinements, including explicit E×B shear suppression and scaled ETG transport, produce quantitative but not qualitative improvements. This work lays the foundation for predictive modeling of future devices. This paper is on arxiv and has been submitted to Nuclear Fusion.« less

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $$\nabla_v \cdot(v\nu(v)f(v))$$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $$\nu(v)$$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational costmore » of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.« less

Microturbulence can produce stationary fine-scale radial corrugations on the plasma density and temperature gradients in magnetic confinement fusion devices. We show that these structures play a significant role in regulating turbulent transport. We focus on the pedestal, studying electron-temperature-gradient (ETG) mode destabilisation and saturation in the presence of radial corrugations on the electron temperature gradient that could result from microtearing turbulence. A linear dispersion relation is derived for a shearless slab case, which indicates that in the presence of a sinusoidal background corrugation, each ETG mode splits into three distinct eigenvalues, with one being the original, one being more unstablemore » and one being less unstable. However, despite the presence of more unstable linear modes, nonlinear gyrokinetic simulations of ETG with corrugated background electron temperature show a reduction of fluxes. Our investigation reveals a radial variation of the phase velocity of the modes that is proportional to the diamagnetic drift velocity and the local pressure gradient. The associated profile shearing breaks the turbulent eddies apart, reducing the transport level. This profile shearing resulting from fine-scale pressure corrugations could be a ubiquitous turbulence saturation mechanism not just in Fusion plasmas, but in Astrophysics and other areas.« less

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking atmore » positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.« less

Extracting reliable information from diagnostic data in tokamaks is critical for understanding, analyzing, and controlling the behavior of fusion plasmas and validating models describing that behavior. Recent interest within the fusion community has focused on the use of principled statistical methods, such as Gaussian process regression (GPR), to attempt to develop sharper, more reliable, and more rigorous tools for examining the complex observed behavior in these systems. While GPR is an enormously powerful tool, there is also the danger of drawing fragile, or inconsistent conclusions from naive GPR fits that are not driven by principled treatments. Here we review themore » fundamental concepts underlying GPR in a way that may be useful for broad-ranging applications in fusion science. We also revisit how GPR is developed for profile fitting in tokamaks. We examine various extensions and targeted modifications applicable to experimental observations in the edge of the DIII-D tokamak. Finally, we discuss best practices for applying GPR to fusion data.« less

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permitsmore » turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Here our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ~40%—through multiscale interactions.« less

Gyrokinetic simulations of turbulence are fundamental to understanding and predicting particle and energy loss in magnetic fusion devices. Previous works have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius effects. This research moves beyond models to demonstrate important corrections using a gyrokinetic Fokker-Planck collision operator with the exact field-particle terms, in realistic simulations of turbulence in magnetically confined fusion plasmas. The exact operator shows significant corrections for temperature-gradient-driven trapped electron mode turbulence and zonal flow damping, and for microtearing modes in a Joint European Torus pedestal under ITER-like wall conditions.more » Analysis of the corrections using parameter scans motivates an accurate model which closely reproduces the exact results while reducing computational demands.« less

One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a novel multi-network physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. Furthermore, this technique presents a novel paradigm for the advancedmore » design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.« less

A reduction of the pedestal pressure with increasing separatrix density over pedestal density (nesep/neped) has been observed in JET. The physics behind this correlation is investigated. The correlation is due to two distinct mechanisms. The increase of nesep/neped till ≈0.4 shifts the pedestal pressure radially outwards, decreasing the peeling-balloning stability and reducing the pressure height. The effect of the position saturates above nesep/neped ≈ 0.4. For higher values, the reduction of the pedestal pressure is ascribed to increased turbulent transport and, likely, to resistive MHD effects. The increase of nesep/neped above ≈0.4 reduces ∇ne/ne, increasing ηe and the pedestal turbulentmore » transport. This reduces the pressure gradient and the pedestal temperature, producing an increase in the pedestal resistivity. The work suggests that the increase in resistivity might destabilize resistive ballooning modes, further reducing the pedestal stability.« less

Recent work on solving partial differential equations (PDEs) with deep neural networks (DNNs) is presented. The paper reviews and extends some of these methods while carefully analyzing a fundamental feature in numerical PDEs and nonlinear analysis: irregular solutions. First, the Sod shock tube solution to the compressible Euler equations is discussed and analyzed. This analysis includes a comparison of a DNN-based approach with conventional finite element and finite volume methods, and demonstrates that the DNN is competitive in terms of degrees of freedom required for a given accuracy. Further, the DNN-based approach is extended to consider performance improvements and simultaneousmore » parameter space exploration. Next, a shock solution to compressible magnetohydrodynamics (MHD) is solved for, and used in a scenario where experimental data is utilized to enhance a PDE system that is a priori insufficient to validate against the observed/experimental data. This is accomplished by enriching the model PDE system with source terms that are then inferred via supervised training with synthetic experimental data. The resulting DNN framework for PDEs enables straightforward system prototyping and natural integration of large data sets (be they synthetic or experimental), all while simultaneously enabling single-pass exploration of an entire parameter space.« less